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ABSTRACT 



Most stars form in compact, dense embedded clusters 
with memberships ranging from a dozen stars to many 
millions of stars. Embedded clusters containing more 
than a few hundred stars also contain O stars that disrupt 
the nebula abruptly. There are observational indications 
that the expulsion of gas from such clusters may occur 
within, or even quicker than, a stellar crossing-time scale. 
Because the star-formation efficiency is empirically de- 
duced to be typically lower than 40 per cent, such clus- 
ters are left severely out of dynamical equilibrium with 
a super-virial velocity dispersion when they are exposed. 
As a consequence, they expand drastically. The velocity 
dispersion of expanding star-burst clusters may achieve 
values of a few tens km/s, which would have a signifi- 
cant impact on the stellar velocity distribution function 
and therefore the morphology of a galaxy. An applica- 
tion of these ideas to the vertical structure of the Milky 
Way (MW) disk suggests that such popping star clusters 
may constitute the missing heat source for understand- 
ing the age-velocity-dispersion relation. If the early MW 
disk experienced a vigorous epoch of star formation then 
perhaps even the thick disk may have formed from pop- 
ping clusters. Thus we have an alternative hypothesis to 
understanding thick disks which does not rely on, or may 
act in addition to, sinking galaxy satellites or sinking cos- 
mological sub-structures. The GAIA mission will test 
this hypothesis by allowing the kinematical fields around 
young clusters to be mapped to exquisite precision. ' 

Key words: Methods: N-body simulations, Astrome- 
try. Stars: formation. Galaxy: kinematics and dynamics. 
Galaxy: evolution, open clusters and associations: gen- 
eral. 



1. INTRODUCTION 



Quoting from the review by Lada & Lada (2003), "obser- 
vations of nearby cloud complexes indicate that embed- 



'To appear in the Proceedings of the Symposium "The Tliree- 
Dimensional Universe with Gaia", 4-7 October 2004. Observatoire de 
Paris-Meudon, France, eds: C. Turon, K.S. O'Flaherty, M.A.C. Perry- 
man (ESA SP-576). 



ded clusters account for a significant (70-90 %) fraction 
of all stars formed in GMCs" (giant molecular clouds). 
This exciting discovery has deep implications for the stel- 
lar distribution function in galaxies, and it is verified at 
the other extreme of the star-formation scale, namely in 
violently interacting gas-rich galaxies, de Grijs et al. 
(2003) find that star cluster formation is a major mode of 
star generation in galaxy interactions, "with > 35 % of 
the active star formation in encounters occurring in star 
clusters". Indeed, the entire Milky Way (MW) popula- 
tion II stellar spheroid could have been build-up from a 
population of star clusters formed during the chaotic early 
merging epoch of the Milky Way (MW), with the present- 
day globular clusters forming the remaining high-mass 
end of the whole distribution (Kroupa & Boily 2002). 

In order to assess the implications that the birth of stars 
in clusters may have on galactic structure and kinemat- 
ics, we need to get a feeling for the physics of star clus- 
ter formation. This, however, leads us rather rapidly to 
the limits of our present means, because the formation of 
star clusters is intimately linked to stellar feedback during 
the first few lO'' yr. This is evident empirically through 
the low star-formation efficiency in cluster-forming cloud 
clumps of less than 40 per cent (Lada & Lada 2003). 
The theoretical problem cannot be solved today because 
radiation transfer in a highly dynamic, turbulent hydro- 
dynamical medium that contains complex chemical net- 
works and magnetic fields is untractable. Todays the- 
ory of star-cluster formation is limited to purely hydro- 
dynamic computations (eg. Mac Low & Klessen 2004; 
Bonnell, Bate & Vine 2003) that give an important in- 
sight into early fragmentation processes but cannot cap- 
ture the essential physics that act to expose the cluster. 
We therefore need to resort to observations to constrain 
the important physical processes responsible for expos- 
ing an embedded cluster. 

Table [2 provides a few well studied examples of very 
young stellar systems, some of which have not yet evac- 
uated their remaining gas. Estimates of the crossing and 
relaxation times are also indicated^, and it can be seen 
that tcr > age for the embedded low-mass objects that 

^The nominal crossing time is ta = 2 Reel /""ccl ~ 
4 (lOOMQ/Meci)^/^ (Heci/pc)^/^ Myr; the relaxation time is 
ticl =0.1 N/lnN tcT- Here Red is the characteristic radius and 
(Tccl is the stellar 3D velocity dispersion; A'^, M^d are the number of 
stars and the stellar mass in the cluster. 
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T-A (emb) 


1 


30 


0.5 


3.6 


4.1 


p Oph (emb) 


1 


100 


0.4 


3.5 


1.6 


ONC (evac) 


1 


5000 


0.5 


18 


0.3 


R136 (evac) 


2 


10^ 


2 


520 


0.6 


Ant. (evac) 


< 10 


10^ 


4 


3600 


0.5 



Table 1. A few well-studied examples of very young clusters: 
Taurus-Auriga (T-A) constitutes perhaps the smallest unit of 
clustered star-formation and is heavily sub-clustered (Bricefio 
et al. 1998). p Oph again is a low-mass cluster that contains 
no massive stars and which is also sub-clustered (Bontemps 
et al. 2001). The Orion Nebula Cluster (ONC) has already 
evacuated its gas and contains no sub-structure (Wilson et al. 
1997; Hillenbrand & Hartmann 1998; Scally & Clarke 2002). 
R136 is the central cluster in the 30 Dor star-forming region in 
the Large Magellanic Cloud (Brandl et al. 1996; Bosch et al. 
2001 ). It already appears to have blown out most of its gas, a 
large part of the cluster having been exposed. It contains no sig- 
nificant sub-clustering. The star-burst clusters in the interacting 
Antennae (Ant.) galaxies have typical radii ^ 4 pc (Whitmore 
et al. 1999). tcr and t,ci ore estimated using R^ci ~ R cind an 
average stellar mass o/0.4 Mq. 



do not contain massive stars, while the opposite is true 
for the rich clusters. The massive clusters are well mixed 
(<cr < age) but not yet relaxed, so some information from 
their formation is still evident. Note that the evacuated 
clusters are likely expanding due to recent gas blow-out 
(Section f2.3t and so R > Red probably. 

The important lesson from this table is that, while the 
radii of the objects are always similar, their masses vary 
by about 5 orders of magnitude. In order to make first 
steps towards understanding the implications star-cluster 
birth may have on galactic scales we therefore need to 
worry about the distribution of Meci; -Rod can be assumed 
constant for now. 



THE BIRTH OF A CLUSTER AND ITS IMPLI- 
CATION 



2.1. The embedded cluster 

Assume a region of a molecular cloud becomes gravita- 
tionally unstable and contracts. Contraction appears to 
occur very rapidly, on a time-scale of 1-3 Myr for a wide 
range of cloud masses and sizes, (Hartmann, Ballesteros- 
Paredes & Bergin 2001). The free-fall collapse time, 
te = for a clump with a radius of 5 pc and a 

mass of 10^(10^) Mq is about 6 (0.6) Myr As the cloud 
clump contracts it fragments because the increasing den- 
sity leads to smaller Jeans masses, Mj cx 1/ ^/p. Turbu- 
lent motions also lead to localised fragments (Mac Low & 
Klessen 2004). Within the fragments proto-stars form on 
a time-scale of 10^ yr (Wuchterl & Tscharnuter2003) and 
decouple from the hydro-dynamical flow thereby becom- 
ing ballistic. They consequently orbit within the contract- 



ing overall cloud clump. Age-dating of the young stars in 
the ONC, for example, suggests that the stars formed over 
a time-span of about 1 Myr (Hillenbrand 1997; Palla & 
Stabler 1999; O'dell 2001; Hartmann 2003), and so most 
of the stars have orbited the forming embedded cluster 
a few times before star formation was halted, probably 
due to the emergence of the O stars. Indeed, the short 
photo-evaporation time-scale of proplyds suggests that 
the central star 91 C Ori with a mass of about 50 
may only be a few 10"* yr old (O'dell 2001). A similar 
picture emerges for the R136 cluster (Massey & Hunter 
1998; Bosch et al. 2001). 

Because most of the stars have orbited a few times un- 
til that critical time when the feedback energy from the 
massive stars ionises the gas throughout the cluster re- 
gion (with a typical radius of a few pc) thereby mostly 
halting further star formation, the embedded cluster must 
be, at that instant, close to dynamical equilibrium. The 
stellar velocity dispersion is 
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where e = Med/ (Med + A/g) sa 20 - 40 per cent is the 
measured star-formation efficiency (Lada & Lada 2003). 
Now note that for 



CToci = 40 km/s, 



M< 



ccl 



eRc 



= lG^-^Mopc-\ (3) 



Since i?cei ~ 3 pc while e w 1/3 we find that an em- 
bedded cluster with the mass of a typical globular cluster 
has a velocity dispersion characteristic of the thick disk of 
the MW, and indeed of dwarf galaxies. The notion thus 
emerges that the birth of massive clusters may have an 
important impact on the structure and internal kinemat- 
ics of galaxies (Kroupa 2002). But this can only be true 
(i) if such large velocity dispersions are actually observed 
in embedded clusters, and (ii) if the embedded velocity 
dispersion can be carried into the galactic field while the 
cluster emerges from its cloud core. 



2.2. Dynamical state just before gas expulsion 

Observations indicate that condition (i) may well be 
fulfilled: 

ONC: age w 1 Myr; a = 4.3 km/s ±0.5 km/s 
(Jones & Walker 1988: proper-motion study, 900 stars) 

Rl 36: age « 2 Myr; cr = 55 km/s ±19 km/s 
(Bosch et al. 2001: spectroscopic l.o.s. velocities, 
48 O and B stars) 

The velocity dispersion for R136 is the three-dimensional 
dispersion, based on the Bosch et al. (2001) line-of-sight 
value CTios = 32±llkm/s. The ONC (Wilson et al. 1997) 
and probably also R136 (Brandl et al. 1996) contain lit- 
tle gas and thus the measured 3D velocity dispersion is 
probably smaller than the pre-gas expulsion stellar veloc- 
ity dispersion within the embedded clusters, a < CTed, 
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since the clusters are probably expanding unless the star- 
formation efficiency was 100 per cent (e — 1). Further- 
more, for the ONC the velocity dispersion is only a lower 
limit because the proper motion analysis removes radial 
and tangential motions. In both clusters the observed stel- 
lar mass is too small to account for the measured velocity 
dispersion. The measured super-virial nature of the ONC 
has been a puzzle for some time (Jones & Walker 1988; 
Hillenbrand & Hartmann 1998). The astonishingly high 
velocity dispersion in R136 is attributed by Bosch et al. 
(2001) to be the result of the orbital motions of massive 
companions. The sample can be split into two mass bins 
each containing 24 stars. Quoting the line-of-sight dis- 
persions and the average mass, 

form > 23.5 M0, m = 49.6M0, cti = 28 km/s, 
for r7i< 23.5 Mq, to = 19.4 M©, cts = 37 km/s. 

If the large a were due to orbital motions then the more 
massive sample ought to have a larger velocity dispersion 
since on average the binaries have more binding energy. 
Instead a\ < (72 which indicates that R 136 is mass seg- 
regated with a kinematically decoupled core of massive 
stars. Mass segregation is indeed evident in star-counts of 
the cluster (Selman et al. 1999), and in a mass-segregated 
cluster the massive stars are expected to have a lower ve- 
locity dispersion. Therefore probably a substantial part of 
a may be due to the motions of stars through the cluster. 

Both clusters thus appear to be well out of equilibrium 
and therefore expanding. Observations of the veloci- 
ties of stars of later spectral type in R136, and improved 
proper-motion studies of the ONC, are needed to further 
test this notion. 



2.3. Gas expulsion 

Concerning condition (ii), CTccI will be conserved to some 
degree i/the star-formation efficiency is small (e < 40 per 
cent) and gas expulsion occurs on a dynamical time or 
shorter (Tg < <cr)- This is not easy to verify empirically 
because we typically do not catch the cluster just at the 
time when it is expelling its gas, although the ^ 0.2 Myr 
old Treasure Chest cluster appears to be doing just that. 
Smith, Stassun & Bally (2004) find the Treasure Chest 
cluster to be very compact with a radius of less than 1 pc. 
The HII region is expanding with a velocity of approxi- 
mately 12 km/s such that the cluster volume will be ex- 
cavated within a few 0.1 Myr. Furthermore, the star-burst 
clusters in the Antennae have t^r ~ 0.5 Myr and Whit- 
more et al. (1999) and Zhang, Fall & Whitmore (2001) 
find gas outflow velocities of 25-30 km/s. This corre- 
sponds to a gas evacuation time-scale of 0.2 Myr which 
is comparable to ici ■ The star-burst clusters in the Anten- 
nae galaxies would thus appear to support explosive gas 
expulsion. Taking the evidence from the ONC and R136 
into account, it appears that in the presence of O stars 
explosive gas expulsion may drive early cluster evolution 
independently of cluster mass. 

Another handle on the evacuation time-scale can be ob- 
tained by considering the energy input from massive stars 
and comparing this to the binding energy, i?bin, of the 
nebula. We have 
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8.6 X lO^s 


0.48 


10^ 


8.6 X IQSO 


0.15 



Table 2. Binding energy and crossing time of two characteristic 
embedded clusters. 

I-Ebini = ° = 8.6 X 10" eig, 

and some characteristic values are listed in Table |2] for 
i?cci = 1 pc. 

Resorting to Maeder's (1990) stellar-evolution tracks 
to evaluate the mechanical (wind) and radiation energy 
output by massive stars we find that a star with mass 

TO = 15 Af© injects 3 x 10^" erg per 0.1 Myr, and 
TO = 85 Mq injects 3 x lO'^^ erg per 0.1 Myr 

The energy injected within a time shorter than the cross- 
ing time is thus larger than the binding energy of the neb- 
ula, and this supports the conjecture that nebula disrup- 
tion may typically occur within a dynamical time-scale, 
i.e very violently, explosively. 



2.4. Evolution of the exposed cluster 

As a result of explosive expulsion of about 70 % of a 
cluster's mass the cluster stars are left orbiting with super- 
virial velocities and expand nearly freely outwards. The 
classical energy argument suggests that such clusters can- 
not survive (Hills 1980; Boily & Ki-oupa 2003). Ba- 
sically, Pleiades- and Hyades-type open clusters ought 
not to exist, unless the stellar IMF with which they were 
formed had no O stars (Elmegreen 1983; Kroupa, Aarseth 
& Hurley 2001, hereinafter KAH). This was a very inter- 
esting proposition, but subsequent observations showed 
that clusters with N > 500 also contain massive stars; ev- 
idence for a truncated IMF has never been found. 

This problem of not being able to understand how open 
clusters form in view of (1) explosive gas expulsion and 
(2) 0.2^ej50.4, has recently been solved by the ad- 
vent of TV— body work which treats the energy exchanges 
between stars due to gravitational encounters precisely. 
Such work has been made possible through theoreti- 
cal research by Sverre Aarseth at Cambridge and Seppo 
Mikkola at Turku on finding numerical and arithmetical 
solutions to the non-linear gravitational dynamics prob- 
lem (Aarseth 1999).^ Using an amended version of the 
Aarseth integrator Nbody6 on a stand-alone PC, KAH 
calculated the evolution of an initially dense embedded 
cluster containing 10'' stars and brown dwarfs, with a 
100 % fraction of binaries. The initial ONC model was 

' Recent work in this field has produced special-purpose gravity- 
pipe-line (Grape6) supercomputers that can achieve a peak perfor- 
mance of 1 TFLOP per sec (Makino et al. 2003), while the next gen- 
eration Grapes machines are expected to reach a peak performance of 
the order of peta-FLOPS. 
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taken to be the most likely configuration from a param- 
eter survey of all possible dynamical states of the ONC 
(Kroupa, Peti- & McCaughrean 1999; Kroupa 2000). Fol- 
lowing the seminal work by Lada, Margulis & Dearborn 
(1984), KAH modelled the gas component (initially 2/3rd 
of the total mass) as a time-evolving background poten- 
tial. This approximation makes the calculations feasi- 
ble on modern computers and is nicely consistent with 
the much more CPU intensive SPH modelling that also 
needs to treat radiative energy transfer in a highly sim- 
plified way (Geyer & Burkert 2001). The KAH models 
include stellar-evolution-induced mass loss and a solar- 
neighbourhood tidal field, and thus constitute the most re- 
alistic existing computations of all relevant physical pro- 
cesses acting during the emergence of a young cluster 
from its natal cloud. 

The overall evolution of the embedded cluster with 
O stars is visualised in Fig. Q which shows the reaction 
of the mass shells to sudden gas removal. The Lagrange 
radii remain constant during the embedded phase, but 
during the final stage just before gas expulsion some mass 
segregation becomes apparent through a contraction of 
the core radius. Gas expulsion leads to expansion of the 
Lagrange radii, and by 1 Myr the theoretical velocity dis- 
persion has dropped to the observed value in the ONC. 
The radial density profile also matches the observed ONC 
by this time. The 50 % mass radius expands with a ve- 
locity of about 1.2 km/s. It crosses the tidal radius by 
about 1 Myr by which time about 50 % of the embedded 
cluster mass becomes unbound. The inner « 25 per cent 
Lagrange radius as well as the core radius contract after 
a few Myr to reach a new minimum by 10 Myr. 

A young open cluster has thus formed. It retains ap- 
proximately fst = l/3rd of the population formed in 
the embedded cluster, and it fills its Roche lobe en- 
tirely. By 100 Myr it evolves to an object that resem- 
bles the Pleiades. Theoretical cluster evolution models by 
Portegies Zwart et al. (2001) furthermore show that the 
Pleiades will evolve to a Hyades-like cluster, and so we 
are left to marvel at the amazing coincidence that within 
less than one-third of a quadrant on the sky we can see 
the following evolutionary sequence with the naked eye: 

ONC =^ Pleiades =^ Hyades 

« 99 Myr w 500 Myr 

Much more work needs yet to be performed in order to 
better understand the role of two-body relaxation and 
binary-star energy exchanges during the critical gas- 
expulsion phase, and to quantify the function 

fst = /st ( Mod , /bin , e , Tg ) , (4) 

where /bin is the primordial binary proportion, and to 
include the presence of the molecular cloud behind the 
ONC which is Ukely to affect the kinematical field. Nev- 
ertheless, these results have essentially solved the open- 
cluster-formation problem. Bound clusters form readily 
despite a low star-formation efficiency and despite gas 
expulsion on a dynamical time-scale. 



Figure 1. The evolution of Lagrange radii (from bottom to top: 
radii containing 5, 10, 20... 50 % of the stellar mass), of the 
gas mass (Mg), of the core radius (Rc), and of the tidal radius 
(Rtid), of an ONC-type cluster. The initial embedded model 
has a central density o/lO'*'* stars/pc'^ with (Joci = 6.8 km/s, 
blows-out its gas after 0.6 Myr within a crossing time, and re- 
sembles the ONC by 1 Myr. After re-virialisation of about l/3rd 
of the initial population the open cluster matches, by 100 Myr, 
the Pleiades (from KAH). 

2.5. Predictions 

These results suggest that each cluster probably forms 
with about three times as many stars than are contained 
within the later re-virialised cluster. There are two mov- 
ing groups (MGs) associated with each cluster, 

MG I: the "new" MG stemming from explosive 

gas expulsion, and 
MG II: the classical MG stemming from cluster 

evaporation, 

such that the relative population of each is given roughly, 
at some time t after the formation time of the embedded 
cluster, by 

Nmgu < Nmgi, for t < 
A^MGii « (l/3)iVMGi/or t> 

where 

Te,/[Myr] « 10(Af/[Mo])°-75 

(Baumgardt & Makino 2003) is the evaporation or disso- 
lution time of a secularly evolving open cluster of mass 
M on an approximately circular Galactic orbit close to 
the solar radius. The velocity dispersion of MGII is 
CT « km/s because the stars leak out of the cluster with 
an energy close to zero (ignoring the rare runaway stars 
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that are shot out as a resuh of three- or four-body en- 
counters near the cluster core). The velocity dispersion 
of MGI is a function of the embedded configuration, 

cmgi = ICth{Mcc\, e,Tg) < (Toci, (6) 

and this can be transformed to a correlation between 
the observed velocity dispersion and the observed clus- 
ter mass, Afci, at some time t <C Tov, 

CTMGI = A^obsl^'^cl), (7) 

which is expected to hold true for an ensemble of MGIs 
and for average values of e, Tg. We might thus be able to 
constrain important star-formation parameters related to 
feedback energy by studying eq. ^theoretically and eq. 
observationally. 

The theoretical research is time-consuming because the 
parameter space spanned by Mcci, i?oci, e, Tg, /bin needs 
to be sampled using precise iV-body calculations. For 
Moci^IO'^Mq (>10^ stars) this work will need the 
afore mentioned special-purpose Grape6 supercomput- 
ers (Baumgardt & Makino 2003). The observational re- 
search will need space-based all-sky astrometry missions 
that capture a substantial fraction of local stars. 



2.6. The role of GAIA 

According to the above theory, each young cluster should 
be surrounded by a radially expanding population of stars 
of similar age as the cluster population. The expansion 
field will not be exactly spherically symmetric because 
variations in the tidal field deflect the stellar orbits. Nev- 
ertheless, a relation of the form of eq. should emerge 
upon measurement of the kinematical field surrounding 
young clusters. 

The GAIA mission is ideally suited to uncover relational 
locally in the MW and in the Large and Small Magel- 
lanic Clouds (respectively LMC, SMC). Locally the high 
precision of the proper motion measurements will allow 
highly detailed mappings of all existing MGI and MGII. 
In the distant LMC only the early-type stars can be mea- 
sured around young clusters, but nevertheless /Cobs ought 
to become evident. 

The tangential velocity to the line-of-sight is 

VT[km/s] = 4.74 x Ai[as/yr] x d[pc], (8) 

where /i is the proper motion in arc sec per year, and d 
is the distance of the star Neglecting pre-main sequence 
evolution for the sake of illustration, an A type star with 
My = LI at the distance of the LMC (d = 50 kpc) 
has an apparent magnitude my — 19.6. From table 8 
in Ferryman et al. (2001) the measurement accuracy for 
such a star corresponds to 28 km/s. This should uncover 
radial flows away from massive young clusters if young 
stellar samples can be combined in sub-regions around 
the clusters. A G-type star with My = 5 at a distance of 
the ONC (450 pc) has my = 13.3. For such a star GAIA 
will have a precision of 0.01 1 km/s, which will suffice to 
define the kinematical state of the entire ONC, although 
issues of crowding do arise. 



3. CREATING A GALAXY 



The ideas developed above can be explored further theo- 
retically by considering their implications on galaxy mor- 
phology. In Section lzTl it was already noted that the birth 
of massive clusters may impact the vertical structure of a 
disk galaxy such as our own MW. In order to study this 
in more detail it is necessary to build a field population 
from an ensemble of popping clusters. 



3.1. Basics 

The basic equation for the velocity distribution function 
of stars is obtained by adding all expanding populations 
from all clusters formed in one "epoch". Note that Wei- 
dner, Kroupa & Larsen (2004) find that the ECMF ap- 
pears to be completely populated for star-forming epochs 
lasting typically about 10 Myr, independently of the star- 
formation rate (SFR). For this purpose we need only to 
consider the distribution of the clusters by mass formed 
in each epoch since the radii of very young clusters do 
not differ significantly (Section^, 

I?(i;z;Mocl,inax,/^) = 

Aficl max 

/ i?(^;,;M,ei)Af(M,ei)ecci(M,ei,/3)dMeci, (9) 

where Vdvz is the resulting number of stars with z- 
components of their velocity vectors in the interval 

V^.Vz + dVz, AfocLmin ~ 10 Mq, and Med, max IS 

the maximum embedded cluster mass in the cluster en- 
semble. The product N^ecidMed is the number of 
stars formed from clusters with masses in the range 
Meci, Mod + dMcci- The expanding population of MGI 
stars emanating from a cluster of mass Mod has a ve- 
locity distribution function D, where D{vz; Afod) dvz is 
the number fraction of cluster stars with velocities in the 
range Vz,Vz + dvz- The physics of star formation is con- 
tained in D{vz; Med), but for an assessment of the im- 
pact of clustered star formation on vertical Galactic disk 
structure it suffices to consider a Gaussian Vz distribu- 
tion, i.e. a Maxwellian speed distribution in the vertical 
direction, with a dispersion given by eq. |2l The actual 
D(w2; Mod,max) would need to account for the MGII, 
i.e. the fraction of stars that re-virialise in the bound clus- 
ter, and which, upon secular evaporation from the cluster, 
contribute a population of stars with a small velocity dis- 
persion to v. Gravitational retardation would also need 
to be incorporated; the stars that expand from a cluster 
are decelerated by the mass within their distance from the 
cluster, such that ctmgi < f ooi- These points are incorpo- 
rated in the exploratory investigation of Kroupa (2002). 
Before proceeding further we need to specify the embed- 
ded cluster mass function. 



3.2. The mass function of embedded clusters 

As is apparent from eq. |9] the fundamental distribution 
function that governs V is the embedded cluster mass 
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function (ECMF). Most empirical constraints suggest the 
ECMF to be a power-law, 

^ecl (X M-f , (10) 

dNcci = Ccci dMcci being the number of embedded clus- 
ters in the mass interval Mod, Afcd + dMcd- 

Lada & Lada (2003) find, for the local embedded cluster 
sample with masses between about 20 and 10'^ Mq, that 
the ECMF can be described as a power-law with /3 k, 
2, with a possible flattening below roughly 50 Mq. The 
form of the ECMF below about a few dozen Mq is very 
uncertain though because such stellar groups are difficult 
to find. They disperse within a time comparable to the 
gas removal time-scale (Adams & Myers 2001; Kroupa 
& Bouvier 2003). 

The maximum mass of a cluster, Afcd.max, in a freshly 
hatched cluster ensemble correlates with the SFR of a 
galaxy (Weidner et al. 2004). In order to extend the 
ECMF to larger masses we therefore need to study galax- 
ies that have a sufficiently high SFR. The Large and Small 
Magellanic Clouds are such systems, and Hunter et al. 
(2003) find 2 < /3 < 2.4 for 10^ < M^^i/Mq < 10*. 

The Antennae galaxies are composed of two currently 
merging major disk galaxies, and the SFR in this sys- 
tem is sufficiently high to sample the cluster MF to 
10'^~^ Mq. For a sample of < 10 Myr old clusters in 
the Antennae galaxies, Zhang & Fall (1999) arrive at 
P = 1.95 ± 0.03 for clusters with masses in the interval 
10^ < Moci/Mq < 10^ and ages 2.5 < t/[Myr] < 6.3. 
"Clusters" more massive than about 10^~ ' Mq are ob- 
served to be composed of star-cluster complexes (Kroupa 
1998; Zhang et al. 2001), thus setting an empirical maxi- 
mum mass of about 10^^^ Mq for a true star cluster, i.e. 
an equal-age, equal-metallicity population (Weidner et al. 
2004). Ki-oupa (1998) and Fellhauer & Kroupa (2002a,b) 
propose such cluster complexes to be the precursors to a 
number of recently found exotic objects: ultra-compact 
dwarf galaxies and low-density clusters. 

For a variety of environments the ECMF thus appears to 
have 

(3^2. (11) 

This is rather remarkable suggesting a possible universal- 
ity of the ECMF 

Note that, apart from the Lada& Lada work, the above 
studies actually report the shape of the mass function 
of young clusters rather than embedded clusters. The 
studied cluster systems typically have ages younger than 
10 Myr Most clusters younger than 10 Myr would still 
contain the majority of their stars within a relatively com- 
pact region even if they are expanding as a result of gas 
blow out (Fig. 0. Nevertheless, this is an issue worth 
keeping in mind; the ECMF underlying the young-cluster 
systems of Hunter et al. (2003) and Zhang & Fafl (1999) 
may well have a somewhat different jS. Kroupa & Boily 
(2002) stress that as a result of the low e a distinction 
between the cluster mass function and the ECMF is nec- 
essary, because clusters loose significant stellar mass as a 
result of gas expulsion. Clearly, this is an important issue 
for future iV-body work. 



3.3. Thickening of galactic disks through clustered 
star formation 

Given the ECMF and a pre-scribed Med, max, eq. |9lcan 
be calculated. For /3 = 1.5 and different Mod, max, the 
resulting T>{vz) is non-Gaussian, having a "cold" peak 
with broad warm to hot wings (Fig.|2ji. The width of the 
wings is determined by Mod. max, while the strength of 
the cold peak is given by Med,min and the cluster-cluster 
velocity dispersion. The distribution of the z-components 




-400 -200 200 

V, [pc/Myr] 



Figure 2. The resulting distribution function of vertical ve- 
locities for a population of stars formed in one epoch of star 
formation which produced a power-law ECMF with j3 = 1.5. 
The distribution function is shown for different maximum clus- 
ter masses, Afcci.max ~ 10*, 10^ "\ 10^ Af0. The curve la- 
belled "G" is a Gaussian distribution function with a velocity 
dispersion equal to the Mod, max = 10*'^ Mq case (26 km/s). 
The cluster-cluster velocity dispersion is 5 km/s. Note that 
1 km/s ~ 1 pc/Myr. From Kroupa (2002). 

of the velocity vectors of solar-neighbourhood M dwarfs 
indeed shows non-Gaussian features (Reid, Hawley & 
Gizis 1995) that are, at least qualitatively, consistent with 
this scenario. Again, the GAIA mission will quantify the 
shape of the velocity distribution of field stars allowing 
tests of the present theory. 

An immediate application of this theory can be seen in 
the hitherto unsolved problem that the observed velocity 
dispersion of solar-neighbourhood stars increases more 
steeply with age than is predicted from secular heating 
of the stellar orbits (Fig. |3}- According to the standard 
picture, stars are formed with a small velocity dispersion 
of about 10 km/s. The velocity dispersion increases with 
time because the stars scatter on giant molecular clouds 
and spiral density waves, and because the disk mass in- 
creases through accretion. Jenkins (1992) calculated this 
diffusion of stellar orbits with time (the solid curve in the 
upper panel of Fig. |3}- The observed vertical velocity 
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(Fuchs et al.'Ol) 

:^Thick disk 
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At [10^ yr] 

Figure 3. Upper panel: The age-velocity-dispersion data for 
velocity components perpendicular to the MW disk for solar 
neighbourhood stars as compiled by Fuchs et al. ( 2001 ) are 
shown as solid dots. Jenkins ' ( 1992 } theoretical relation is plot- 
ted as the solid curve. The thick-disk datum is indicated by the 
old star Lower panel: The maximum cluster mass, Meci,max. 
in the population of star clusters needed such that the differ- 
ence in the top panel between observed data and theory van- 
ishes is calculated and plotted as a function of time assuming 
P — 1.9. The shaded region comprises the uncertainty range on 
Afcci.max given by the observational uncertainties on (upper 
panel). The region between the old horizontal dashed lines indi- 
cates the Afeci.max needed to account for the vertical thick-disk 
kinematics, the thick disk having been born between about 10 
and 12 Gyr ago. From Kroupa (2002). 



dispersion data (Fuchs et al. 2001) are plotted in the top 
panel of Fig.|3] The more recent work of Nordstrom et al. 
(2004) shows similar results. It can be seen from Fig. 13 
that the theory cannot account for the observed increase 
of cr^. A "heat source" appears to be missing. 

Adding popping star clusters may resolve this problem 
of the missing heat source. The theoretical velocity dis- 
persion, (Jz, can be calculated readily from I?(fz) as a 
function of M, 



cci.max 



for various (3 and an assumed small 
cluster-cluster velocity dispersion of 5 km/s. As a first 
ansatz we thus say that the clusters form in a very thin 
MW disk. Fig.|3lindicates what sort of Moci.max(i) with 
P — 1.9 would be needed to explain the observed veloc- 
ity dispersion assuming Jenkins' (1992) model applies. 
The result is that Afoci,max takes reasonable values and 
that it decreases towards younger ages. 

The interesting suggestion coming from this approach is 



that the MW disk may have been quietening down over 
time as its gas reservoir was depleted. It may have be- 
gun forming with a high SFR reaching to Afcd.max ~ 
10^ Mq thereby producing the thick disk about 10- 
12 Gyr ago. Such a high SFR in the early primarily 
gaseous MW disk may have been induced by a tidal per- 
turbation from a nearby neighbour galaxy (Kroupa, Theis 
& Boily 2004). Or it may be the result of an early disk 
instability (Noguchi 1999). A full solution to the ECMF 
needed to match the vertical thick-disk kinematics can be 
found in Kroupa (2002). The thin disk may then have 
begun its assembly from star-cluster populations with a 
maximum cluster mass decreasing to the present time. 
The decreasing Afcd.max would come about from a de- 
creasing SFR (Weidner et al. 2004) as the conversion of 
gas to stars proceeded. 

The above constitutes a simple example of how clustered 
star-formation may affect the morphology of galaxies. 
The results are very encouraging, in that the constrained 
parameters Afcci.max, /? are reasonable. Combined with 
the observational evidence discussed in Section|2lthe no- 
tion that star clusters are fundamental building blocks of 
galaxies is sharpened. Clearly, we need to calculate the 
velocity distribution function of expanding stellar popula- 
tions, VlvxjVy, Vz), within a realistic and time-evolving 
Galactic potential in order to construct three-dimensional 
velocity ellipsoids at different locations in the MW disk 
and at some given epoch for a pre-scribed star-formation 
history. This needs treatment of the diffusion of orbits 
in the time-evolving potential. Furthermore, ancient thin- 
disk clusters with Al, 



ccLmax 



<; 10^ Mq may be evident in 
the MW disk today as evolved cluster remnants (eq.|5}. 



4. CONCLUSIONS 



The measured super-virial velocity dispersion of ONC 
and R136 stars suggests that both clusters may be out- 
of-dynamical equilibrium after blowing out their residual 
gas explosively. Explosive gas expulsion from massive 
young clusters is also supported by measurements of the 
outflow velocities of gas, and the Treasure Chest clus- 
ter may be caught in the act of doing so. Young clusters 
containing O stars therefore may expand rapidly after gas 
blow-out loosing a large fraction of their stars. This ex- 
panding population forms a moving group (MGI) with a 
velocity dispersion, cmgi^ that is related to the pre-gas 
expulsion velocity dispersion of stars in the embedded 
clusters, and thus to the mass, Afd, of the re-virialised 
cluster that re-contracts as the nucleus of the expand- 
ing MGI. GAIA will measure this correlation (Jmgi = 
/Cobs(A/ci), through which we will get an important han- 
dle on the basic physics acting during cluster formation. 

Since most stars form in embedded clusters it follows 
that cluster birth may leave non-negligible imprints in the 
kinematical and structural properties of galaxies. An ex- 
ample has been presented by calculating the vertical ve- 
locity dispersion of a population of stars born in an en- 
semble of star clusters. Relating this dispersion to the 
missing heat source in the age-velocity-dispersion rela- 
tion of solar neighbourhood stars allows us to estimate 
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the evolution of the embedded-cluster mass function with 
time. The result of this exploratory work is that the MW 
may have been quietening down somewhat, that early 
epochs of the MW disk may have been characterised by 
more violent star formation within a thin disk producing 
clusters reaching to masses similar to those of present- 
day globular clusters. A thick disk may thus be a natural 
outcome of such events. The merging of cold-dark matter 
sub-halos may not be needed to explain the existence of a 
thick disk, but satellite mergers are of course not ruled out 
as having shaped the old MW disk (Abadi et al. 2003). 

A more detailed theoretical study of these processes is 
needed to quantify the stellar velocity ellipsoid in depen- 
dence of various evolutionary models of the MW disk, 
and to construct statistical models of the MW for future 
comparison with GAIA data. 

Because this notion that popping clusters may shape en- 
tire galaxies is so important, it will also be urgent to ob- 
serve the velocity dispersion and the gas content in as 
many very young clusters as is possible. While GAIA 
will have the defining impact on these ideas, it will not 
be able to measure the kinematics of very young stars in 
dense clusters, so that ground-based proper-motion and 
line-of-sight velocity observations are definitely needed 
to quantify the dynamical state of as large a number of 
very young clusters as is possible. 

The beauty of this entire approach lies in the realisation 
that small-scale baryonic processes may affect galaxy- 
scale structures through stellar dynamics. Star-formation 
scales may thus impact cosmological issues in a hitherto 
unappreciated manner. 
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